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Abstract 

We study a model for the evolutionarily stable strategy (ESS) used by biological populations 
for choosing the time of life-history events, such as migration and breeding. In our model we 
accounted for both intra-species competition (early individuals have a competitive advantage) 
and a disturbance which strikes at a random time, killing a fraction 1 — p of the population. 
Disturbances include spells of bad weather, such as freezing or heavily raining days. It has 
been shown in [23] . that when p = 0, then the ESS is a mixed strategy, where individuals 
wait for a certain time and afterwards start arriving (or breeding) every day. We remove the 
constraint p = 0 and show that if 0 < p < 1 then the ESS still implies a mixed choice of 
times, but strong competition may lead to a massive arrival at the earliest time possible of a 
fraction of the population, while the rest will arrive throughout the whole period during which 
the disturbance may occur. More precisely, given p , there is a threshold for the competition 
parameter a, above which massive arrivals occur and below which there is a behaviour as in 
[23] . We study the behaviour of the ESS and of the average fitness of the population, depending 
on the parameters involved. We also discuss how the population may be affected by climate 
change, in two respects: first, how the ESS should change under the new climate and whether 
this change implies an increase of the average fitness; second, which is the impact of the new 
climate on a population that still follows the old strategy. We show that, at least under some 
conditions, extreme weather events imply a temporary decrease of the average fitness (thus an 
increasing mortality). If the population adapts to the new climate, the survivors may have a 
larger fitness. 

*The first two authors contributed equally to this work. 
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1 Introduction 

Proper timing of life-history events, like emergence, germination, migration or breeding, is crucial 
for survival and successful reproduction of almost all organisms. Timing may be set by endogenous 
rhythms or by extrinsic environmental clues (e.g. day length or temperature, m), but in almost 
all cases timing seems to have evolved according to two contrasting selective pressures. On the 
one hand, the first individuals that emerge or arrive at a given site often perform better, because 
they can profit from the better habitats and benefit from reduced competitions (at least for some 
time). On the other hand, however, early individuals may suffer from higher mortality, as they 
expose themselves to the risk of adverse environmental conditions, which usually are more likely 
early than late in the season. Autumn migration may seem an exception to this pattern, as the 
risk of mortality is probably larger for late than early departing individuals. However, in this case, 
migrants may benefit from a longer stay in their breeding grounds (allowing to rise a further brood 
or acquire larger fat reserves for migration). At the end, this pattern can be seen as the exact 
reverse of the process going on in spring, and can therefore be modelled in the same way. 

In a seminal work Iwasa and Levin [23] have provided a first theoretical description of how the 
risk of incurring in adverse environmental conditions may shape the timing of life history events of 
a population, as a result of evolution over many generations. Under the assumption that adverse 
conditions (disturbance according to their definition, which we will follow hereafter) strike at a 
random time and are so strong that no individual incurring in the disturbance can survive, the 
authors show that (in many cases) the evolutionarily stable strategy (ESS from here on) is an 
asynchronous choice of times in the population. This asynchronicity has been observed in various 
settings and modelled by several authors (see E], DU, BE], EH, [391 just to mention a few). In this 
paper we extend Iwasa and Levin’s model to the broader scenario where the disturbance is soft, 
meaning that it kills each individual with probability 1 — p (in many cases some individuals in the 
population survive even to dramatic adverse conditions). The model of Iwasa and Levin can be 
seen as a particular case of ours when p = 0. 

Before going into the details of our model, we have to mention that it focuses on the long time 
behaviour of a very large population. Indeed it is implicitly assumed that when the existence of 
mutants with better fitness is theoretically possible, then such mutants will appear and spread across 
the population. This does not take into account the disappearance (in finite populations) of certain 
alleles by mere random factors, a phenomenon which can be studied by means of mathematical 
population genetics (see for instance [15]). Long time behaviour of finite populations can also 
be studied through spatial models, namely interacting particles systems. Space not only adds 
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complexity m , but may also be interpreted as “type”, that is the location of one or more individuals 
can be seen as representing their genotype. For the simplest among this models, the branching 
random walk, much has been done: for instance in EQ El m m ESI S3 one finds characterization 
of the persistence/disappearance of genotypes (seen as locations for the model), on general space 
structures; the same can be found, for some random graphs, in El Eg. Stochastic modelling and 
interacting particle systems have been successfully applied to biology and ecology (see kiii sum 
H3l 1231 [251 [261 ESI UHl S8] just to mention a few). Although stochastic modelling is very interesting 
and complex, here we will assume that over many generations, our populations have been sufficiently 
large to justify the use of a model where stochasticity appears only in the random time at which 
the disturbance strikes. 

As for terminology, in the present work the life-history events under study are arrival times, 
meaning that we focus on migratory birds and the time of arrival to their breeding grounds. This 
choice of words should not be considered a reduction in the scope of this paper, since our modelling 
approach is very broad, as it applies to the investigation of the timing of any life-history event when 
the benefits from being early and the risk of incurring in a soft disturbance are in conflict. Migratory 
birds are a well studied biological system where timing is crucial for the fitness of individuals, and 
where a long record of adverse conditions killing or impairing the reproduction of individuals exists 
(see e.g. [33]). Moreover the interest in the timing of these recurrent biological events is very strong 
among biologists after that several studies have consistently observed an advancement in arrival 
and reproduction of birds supposedly as a consequence of climate change (see m, no). 
Climate change not only implies warming temperatures, but also higher frequency and intensity 
of extreme meteorological events (see EH). This increased weather unpredictability may severely 
affect migrant birds, because warmer springs prompt birds towards earlier arrivals, while more 
frequent unseasonable weather increases the risk of mass mortality events (see [2] for a stochastic 
model for random catastrophes striking a spatially structured population). It is widely accepted 
that climate change is endangering migrant populations (see m)- Ornithologists therefore strongly 
need models investigating the contrasting forces affecting the timing of bird migration to improve 
their understating of the ongoing ecological processes and to plan better conservation strategies for 
declining migrant populations. We point out that even if our primary interest here are evolutionarily 
stable strategies which arise in large populations after many generations of stable climate (an 
equilibrium situation), our study also allows us to analyze some effects of a sudden climate change 
(an off-equilibrium dynamics). Indeed our results not only describe the ESS (Theorem 3.5) but 
also the effects of the climate change on the fitness of a population not yet adapted to the change 


(Propositions 3.7, 3.8 and 3.9). 

Here is a short outline of the results of the paper. In Section [2] we introduce the model and 
the notation. We consider an intra-species competition regulated by a parameter a > 0 and we 
define the fitness as a function of the arrival time conditioned on the disturbance time, the 
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expected fitness (/>„ depending only on the arrival time (averaged over all admissible disturbance 
times) and, later on, the average fitness A M (averaged over all disturbance and arrival times). We 
describe the meaning of an evolutionarily stable strategy (ESS) and we discuss the easiest cases 
where either there is no competition (a = 0 ) or the probability p of surviving the disturbance is 1 , 
i.e. the disturbance has no effect whatsoever (Remark |2.1[ ). Section [3] contains our main results. 
Theorem 3.5 extends [23], Appendix A] and shows that there is only one possible ESS for the 


population in response to a fixed disturbance distribution, which we imagine supported in \tf,tf]. 
The behaviour of the ESS depends on the value of ( a,p ). Indeed there exists a function of p, say 
aM = clm(p ), such that if a < aM then the population starts arriving after a date x c > tf and 
there are arrivals every day until tf. if a = om then continuous arrivals starts at tf. if a > aM 
then a fraction 7 of the population arrives at time 0 and the rest arrives continuously starting 
from tf. The dependence on a and p of the fitness of individuals following the ESS is discussed; in 
particular we show (see Remark 3.6) that the theoretical maximum value of the average fitness A^ 
is not attained by any ESS. This proves that, even though an ESS is a strategy that each member 
of the colony considers fair, it is not the best choice for the colony as a whole. The dependences of 
every relevant coefficient on a, p and the disturbance distribution are summarized in a table before 


the beginning of Section 3.1 As an example, we describe the case where the disturbance strikes 


according to a uniform distribution (see Section 3.1). The effects of climate changes are studied in 


Section 3.2 The main questions discussed here are the following. How does the ESS change after a 
climate modification? What happens if a population would keep the same strategy after a climate 
change? A worse climate, that is a smaller p, may change the shape of the ESS delaying the first 
arrivals and increasing the fitness of each individual (provided that the population survives the 
transition). If the distribution of the disturbance is linearly rescaled, the same happens to the ESS. 


The behaviour of the fitness, before the strategy adapts, is studied in Propositions |3.7[ 3.8 and 


3.9 A consistent delay of the disturbance reduces the fitness of all individuals following the former 


ESS. If the disturbance arrives earlier than in the past, then the average fitness of the population 
increases. When the competition is weak then a decrease of p implies a killing of a larger fraction 
of the population and hence a lower average fitness. In Section [4] we discuss and summarize the 
conclusion of the paper. Section [5] is devoted to the proofs of our results. 


2 The model 

Iwasa and Levine |23| studied different ways to model the fact that individuals choosing an early 
date of arrival, if no disturbance were present, would obtain a larger fitness than those arriving 
later. This may be due to the fact that a decrease of reproductive success (better resources at 
earlier times, |23L Case 1]) or to competition between individuals, for instance those who arrive 
earlier may feast on food, while those arriving later will not ( |231 Case 2 ]). We focus on the second 
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case, where Iwasa and Levin proved the emergence of a mixed strategy for the arrival dates (which 
is more interesting and realistic than Case 1, where all individuals choose the same date). 

Suppose that p is the probability measure, supported on [0, +oo), according to which individuals 
choose the date of arrival to the breeding sites. Its cumulative distribution function F fl is defined, 
as usual, by F^{x) := /r((—oo,x]) for all i£l. The population is struck by a single disturbance 
(e.g. storm, frost) whose date is randomly distributed on (ty,f/) with density /. It may be that 
tf = Q(tf is the first possible date of disturbance while tf is the last). Individuals already present 
when the disturbance occurs, survive with probability p 6 [0,1]. The fitness of an individual 
is a decreasing function of the fraction of individuals (in the whole population) who are already 
present when it arrives, and of the parameter a > 0 which represents the strength of intra-species 
competition. Inspired by |23j we choose the following expression for the fitness of an individual 
arriving at date y, given that the disturbance strikes at time x: 


Vv(yk) 


p exp (-aF^y)) if 0 < y < x; 

exp (-apF^x) - aF^y) + aF^x)) if y > x. 


Thus ipn(y\x) is the fitness (under the strategy p) of an individual arrived at time y when the 
disturbance strikes at time x. It is worth noting that we are implicitly exploiting the Law of Large 
Numbers (LLN), since the exact value of ip^{y\x), for instance if y < x, is 


\ x ) = pexp(-aN(y)/N), 


where N(y) is the random number of individuals arrived before y and N is the population size. If 
N is large, the LLN implies that N(y)/N can be approximated by F lt (y). Similarly one proceeds 
with the case y > x. 

We consider the expectation of the fitness, with respect to the disturbance date: the expected 
fitness of an individual arrived at y is 


rtf 

<f>(i(y) = / iptj,(y\x)f(x)dx. 

Jo 

More explicitly, for y > 0, 


ry 


f 


<l>n{y) = exp(-aF M (y)) / exp(a(l - p)F li (x))f(x)dx + p / f(x)dx 

1 Jo Jy 

r f y 

= exp (~aF^(y)) / (exp(a(l - p)F fl (x)) - p)f(x)dx + p 

1 Jo 


( 2 . 1 ) 


We note that, in the last integral of equation (2.1) (as well in every integral of a function of type 
k(x)f{x) in the sequel), one could use +oo instead of tf and the value of the integral would be the 
same (as the definition of 0^); moreover if y > tf, by we mean — f y . 
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We assume that, the population follows evolutionarily stable strategies, that is, distributions of 
arrival times which grant no advantage to any particular choice of arrival date. More precisely, an 
evolutionarily stable strategy (ESS) y, is such that no mutant can have an advantage, namely for 
all y E supp(/x) and for all z E M + , (ft nil/) > 4>u( z )- This implies that = A for all y E snpp(y), 

where A := sup{0 At (y): y E M + } < 1. We recall that supp(//) is the (closed) set of y such that 
y(y — e, y + e) > 0 for all e > 0. We will also need to define the essential support Esupp(/) of 
a real function /, which is the support of the associated measure A i->- f A f(x)dx (for instance, if 
/ is continuous then Esupp(/) = {f 0}). We recall that in [23] the probability distribution of 
the disturbance was considered as absolutely continuous with a single peak density / which was 
taken, on [0 ,tf] as a polynomial vanishing at the extrema of the interval. We only assume that / 
is a probability density, supported in [0,i/]. From now on, without loss of generality, we assume 
tf : = maxEsupp(/) and tf := minEsupp(/). 

We are interested in studying y as a function of a and p (and of /, but here / is thought as 
fixed). The extremal cases where either a = 0 or p = 1 are easy to describe. 


Remark 2.1. 1. If a = 0 and p = 1, there is no competition and the disturbance has no effect. 

Then 4>n(y) = 1 for all y > 0; moreover, every y is an ESS (indeed the disturbance has no 
chance to shape an evolutionary response of the species). 

2. If a = 0 and p < 1, there is no competition. From equation 

p) f* f f(z)dz which is non decreasing and continuous and fi^y) = 1 for all y >tf. We see 
here that there is no dependence on y. A probability measure y is thus an ESS if and only 
if supp(y) C [tf, oo). This means that all idividuals will arrive after the last possible date of 
disturbance. 


(2.1) we have, (fJy) = 1 - (1 - 


3. Ifa> 0 andp = l, there is competition and the disturbance has no effect. From equation (2.1) 


we have 4>n(y) = exp(— aF^fy)) which is right-continuous and nonincreasing. Using the same 
arguments as in the proof of Theorem |ff.5| it is straightforward to prove that there is a unique 
ESS, namely y = 5 q. Everybody arrives at the first possible arrival date (indeed there is no 
risk in doing so). 

4- If a > 0 and p = 0, the case has been studied in l23\. Case 2] and can be retrieved as a 


particular case of Theorem 3.5. There exists a critical date x c (depending on a and f) after 
which individuals start arriving according to an absolutely continuous measure y such that 
supp(y) = ( x c ,tf). The expected fitness of each individual is A = 


The interesting case is when a > 0 and p < 1, that is, competition in the population and 
effective disturbance. These are the constraints which we assume thereafter. 
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3 Main result 


We are able to prove (Theorem 3.5) that given p < 1, there exists a critical am, depending only on 
p (not on /) such that: 


1. if a < cim then the ESS is as follows: there is a critical date x c > tf (depending on o, p and /) 
after which individuals arrive continuously while disturbances are possible (see for instance 
Figure [5]). This extends the previously known result for p = 0 since om{ 0) = +00 (see [23i, 
Case 2 ]); 

2 . if a = a,M then the ESS is as before, with x c = tj (individuals arrive throughout the whole 
period of possible disturbance, see for instance Figure [ 6 ]); 

3. if a > clm then the ESS is such that a fraction 7 of individuals arrive at 0, and the remaining 
arrive continuously during the whole period of possible disturbance (see for instance Figure[7]). 

Before stating our main result, we define the quantities clm, x c and 7 . 


Definition 3.1. If p = 0, then clm(p ) := + 00 ; if p £ (0,1), then clm(jp) is the solution to the 
equation 


r exp(ci M ) dz 


yl—p — 


p 


1 

V 


(3.2) 


Note that the solution to equation (3.2) exists and is unique since the l.h.s. is a continuous, 
strictly increasing function of om which vanishes at om = 0 and goes to infinity as om ~> 00 . The 
inverse function of a« will be denoted by pm{cl) (by definition pm{+ 00 ) := 0). Clearly a < om{p) 
(resp. a > om{p)) if and only if p < pm(cl) (resp. p > pm(cl))- 


Definition 3.2. If a > om{p) let x c = x c (a,p,f) := 0. If a < om{p) define x c as the maximal 
solution to the equation 


J f(x)dx = (l -p + 


exp (a) 


d^ 


T-P — 


P 


-l\ -1 


(3.3) 


Note that the solution to equation (3.3) exists and is unique, since the r.h.s. is positive and 
strictly smaller than 1 , while the l.h.s. is a continuous, nonincreasing function of x c which takes 
values 1 and 0 at x c = 0 and x c = tf respectively. By definition if a < om(p), x c £ Esupp(/). 


Definition 3.3. If a < om(p) let 7 = 7 (a,p) := 0. If a > om{p) let 7 be the unique solution to the 
equation 


exp(— 07 ) 


rexp(a) ^ 


exp(a7) Z P p p 


(3.4) 


7 








The solution exists and is unique since the l.h.s. is a continuous, strictly decreasing function 
of 7 which takes values in (1/p, +oo) when 7 = 0 and is equal to 0 when 7 = 1 . Observe that 
7 ( 0 ,0 ) = 0 for all a > 0 . 

Note that, in accordance with [23], if p = 0 then om = + 00 , 7 = 0, and f* f f(x)dx = a/(l + o). 
Let us discuss some general properties of cim, x c and 7 . For details on the proofs, see Section [5} 

Remark 3.4. 1. The map p i->- aM(p) is continuous and strictly decreasing; see Figure [i| By 

using elementary techniques of implicitly defined functions it is not difficult to prove that 
lim p _ > o+ a M (p) = + 00 , lim^- a M (p) = 0. 

2. Given a fixed f, we have that r c and a 1 —> x c are strictly decreasing and left continuous 
everywhere. The map p 1 - 7 - x c is right continuous at po only in the following cases: (a) 
a > a M (po); (b) a < a M (po) and ( x c (a,p 0 ,f ) - e, x c (a,p 0 , /)] C Esupp(/) for some e > 0; 
(c) a = om(po) and tj = 0. Similarly ae> x c is right continuous at a 0 only in the cases: (a) 
ao > a M (p); (b) a 0 < a M (p) and ( x c (a 0 ,p,f) - e, x c (a 0 ,p, /)] C Esupp(/) for some e > 0; 
(c) ao = om(p ) and tj = 0. Again using standard techniques of implicitly defined functions 
we have 


lim x c (a,p,f)<x c (a,0,f)<t f = lim x c (a,p,f)- 
p-s>0+ a—>0+ 

lim x c (a,p, f) = t f = lim x c (a,p,f). 

P^Pm ( a )~ a —>a M (p) 


(3.5) 


3. The function (a,p) 1— > 7( a,p ) is continuous on (0,+00) x [0,1) and p e-)- 7 (a,p) is strictly 
increasing (for all fixed a G (0, +00) ). Moreover r y(a,p ) < p for allp E (0,1) and for all a > 0; 


a i-T- 7 (a,p) is strictly increasing in [aM(p), +00) for all fixed p E (0,1) (see Remark 5.3). See 
Figures and^for some plots of p 1 —> 7 (a,p) and a 1 - 7 - 7 (p,a). Moreover, we have 


lim 7 (a,p) = 0 , 

P-S-0+ 


lim 7 (a,p) = 1 , 

p— >i~ 


lim 7 (a,p) = p, 

a—>00 


where the last limit holds for all p > 0 . 
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Figure 1: p om(p)- 




- p=0.1 

-p=0.2 

- p=0.3 

. p=0.5 


Figure 2: p H > 7 (p, a); a = 3, 2,1, 0.5. 


Figure 3: a 1 —>• 7 (p, a); p = 0.5,0.3,0.2, 0.1. 


Theorem 3.5. Let clm, x c and 7 be as previously defined. Given a > 0, p £ [0,1) and f a 
probability density supported in [0,t/]) (where by definition tf = maxEsupp(/) ), there exists a 
unique ESS p. In particular, p = 7^0 + (1 — 7)17 where u is an absolutely continuous probability 
measure with supp(z') = [x c ,tf] nEsupp(/). The cumulative distribution function F u is implicitly 
defined, on x > x c , by 


r exp(a(l-7)i7(a;)+a7) ^ 


J exp (cry) 

The value of A (the supremum of cfij is 


yl~P _ 


p 


1 f x 

jJ f(y)dy. 


A = A (a,p) = < 


1 


if a < qm 


1+a -p) fr l “’ Tfej 

pexp(— 07 ) if a> om- 


(3.6) 


(3.7) 
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Finally, om is a phase transition value for the competition in the sense that 


a < om => 7 = 0 , x c > tf] 

a — aj\4 = t > 7 — 0 , x c — if ] 

a > om => 7 E (0, (1 — clm/o) A p), x c = 0. 


It is worth noting that, in equation 


we can write F^fx) instead of 7 + (1 — 7 )F v (x). 


Several features of the ESS p can be read from Theorem 3.5. Beside arrivals at 0 (possible only 
when competition is sufficiently strong), individuals arrive only at possible disturbance dates. This 
means that if it is certain that during [s,t] C [tj.tf], no disturbance is possible (i.e. f(x) = 0 for 
all x E [s,i]), then the probability of arrivals in such interval is zero. This is due to competition, 
which advantages early birds: there is no point in choosing to arrive in between s and t, since there 
is no risk in choosing s instead. 

In most real cases, one may assume that / > 0 on [0,t/]. In that case individuals will arrive 
either (1) avoiding the first part of the possible time period (weak competition, see Figure |Tg|) or 


(2) during the whole interval [0,tj] without massive arrivals at 0 (critical competition, Figure 15) 
or (3) during the whole interval with massive arrivals at 0 (supercritical competition, Figure |~i~i]). 

The first arrival time x c is strictly decreasing with respect to p and to a (Remark |3.4[ ); this 
means that strong competition and/or high probability of surviving the disturbance, push arrivals 
to 0. If p is fixed, competition needs to be above the threshold aM(p), in order to have arrivals at 
0. On the other hand, if competition is fixed, only weak disturbances lead to early arrivals. 

As for 7 , the fraction of arrivals at 0, we know that it increases with p. If p > 0, from 
equation (3.6), using 7 (a,p) | P as a -> 00, it follows that F v (x) — > 1 for all x E (. tf,tf ] as 
a — > 00. This means that, as the competition increases, given that an individual does not arrive 
at time 0 (this probability converges monotonically to 1 — p from above) then the probability of 
arriving after x > tf goes to 0, that is, in the limit as a —> 00 the arrival distribution converges 
to p5o + (1 — p)5t (that is, do if tf = 0). Note that this does not happen when p = 0: in that 
case F v {x) = // f(y)dy which implies that F u converges to the cumulative distribution of the 

disturbance arrival time. 

Associated to a given strategy p, there is the average fitness A M := / 0 + °° 4> fl (y)p(dy). In partic¬ 
ular when p is an ESS, then ep^fy) = X for all y E supp(/i) (where A = A (a,p), see ( |3.7| >), hence 


A m = A (a,p). If a < om we can relate A (a,p) to x c : by equations (3.3) and (3.7) we have 


L 


i ’ mdy =i _h^ A(0 , rt p ) 

1 -p J 1 


dz 


A-p _ 


P 


(3.8) 
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It is not difficult to prove that (a,p) i->- A (a,p) is continuous in (0, + 00 ) x [0,1); moreover a ^ X(a,p) 
and p i-)- A (a,p) are strictly decreasing functions (see Remark |5.3[) such that 


lim A (a,p) = 1, 

a— >0+ 


lim A (a,p) 
0 + 


1 

1 + a’ 


lim X(a,p) = 0; 

a—>•+ 00 

lim X(a,p) = exp(—a), 

p— >i - 


hence A(a, 0) = (1 + a) -1 > A (a,p) > exp(—a) = A(a, 1) for all p E [0,1]; moreover A (clm{p),p) = p. 

Thus, when p is an ESS the average fitness A M = A (a,p) is decreasing with respect to p (when 
a > 0), hence, if we think of the fitness as a probability of survival, the average rate of survivors 
is decreasing if the chance of surviving the catastrophe p is increasing. From a biological point of 
view, this model suggests that, in the presence of competition, if the disturbance is weaker (that 
is, p increases) then the ESS pushes the colony towards an early arrival on the site, increasing the 
negative effects of the competition on the fitness (which overcome the positive effects of the weaker 
disturbance). Hence the stronger the disturbance the higher the average fitness corresponding to 
the ESS (in some sense, if we think of the average fitness of the population as its “strength”, a 
strong disturbance will select a stronger population). 


a-»A(a,0.1) 

A 



0.0 0.2 0.4 0.6 0.8 1.0 P 


a->A(a,0.3) 

A 



0.0 0.2 0.4 0.6 0.8 1.0 P 


a-»A(a.0.8) 

A 



0.0 0 2 0.4 0.6 0.8 1.0 P 


One may wonder if there are strategies p which (given the environment (a,p, /)) lead to a larger 
A m , and whether these strategies are Evolutionary Stable Strategies. 

Remark 3.6. The supremum of the map p (->• f 4> fl (y)p(dy) = X^ is (1 — exp(— a))/a and is 
attained by p if and only if infsupp(/i) > supEsupp(/) = tf and F ^ is continuous (and this 
supremum holds for any fixed p, see details in Section [5]). This means that the best choice for the 
population as a whole, is to start arriving after the last possible date of disturbance. In particular no 
ESS can attain this maximum value (an ESS is supported in [0 ,tf)). An ESS is a strategy which is 
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well accepted by every individual (since the fitness is constant and maximal), while the maximizing 
strategies imply that some individuals “accept” a lower fitness for the benefit of the whole colony 
(note that the continuity of implies that the population cannot choose to arrive simultaneously 
attf, which would guarantee the same fitness for everyone). Moreover, even if the maximal fitness 
is attained only arriving after tf, it is possible to prove that we can find a sequence of strategies 
(not ESSs) {p n } n >i, supported in [0,tj] such that A^„ > (1 — exp(— a))/a — 1/n. In Figure^] the 
dashed line represents sup A , A whereas the dotted and the solid ones are the average fitness when 
the population follows the ESS, A(a, 1) and A(a, 0), respectively. The filled region represents all 
possible values for A (a,p) with p £ (0,1). Note that following the ESS a population cannot achieve 
the maximal average fitness, nevertheless if a is either small or large, then the dashed and the solid 
line are close, hence, if p = 0 or at least p is small, then the average fitness is not so far from its 
theoretical supremum. 



Figure 4: Maximum average fitness and the admissible region for A (a,p). 

In the following table, we summarize the main properties of the coefficients an, 7 , A and x c \ 
by /* x and \ iwe mean that a particular coefficient is increasing or decreasing with respect to 
the parameter x. 
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Coefficients 

Dependence 

Properties 

Phase transition competition clm 

P 

\ P 

CbM G [0, + 00 ] 

Probability of emigration at time 0, -y 

a , P 

/’a /p 

7 = 0 if a < um 

7 G (0 ,p) if a> a M 

Fitness of individuals for the ESS, A 

a , P 

\ a \ p 

exp(—a) < A (a, ■) < (1 + a)^ 1 

Maximum average fitness A for a generic strategy 

a 

A (a,p) = (1 — exp(—a))/a 

Earliest arrival time for an ESS x c 

a, P, f 

\ a \ p 
x c G [tf,tf] if a < aM 
x c = 0 if a > ajf 


Given ( a,p,f ), Theorem |3.5| gives the unique ESS p and its average fitness A. We note that 
the map ( a,p,f ) i —> (A,/x) is not injective. Indeed, at least when a < um(p), it might be that 
A (a,p) = A (b,q) for some q and b. Even if we fix ( a,p ), at least in the subcritical case, different 
disturbance distributions may lead to the same ESS. Indeed, suppose that a < clm{p ) then x c > 0 
and 7 = 0 . Given p is fixed, from equation (3.6), / is uniquely determined on [x c ,tf\. Nevertheless, 
on [0,x c ] the only constraint is f^ c f(x)dx = (p — A)/(l — p), which can be satisfied by infinitely 
many distributions. This means that different levels of competition, and/or of climate, can lead to 
the same response p and same fitness A. 


3.1 Uniformly distributed disturbances 

In this example we suppose that /(f) = (tf — that is, the law of the disturbance is 

uniformly distributed in the interval [tf, tf] (we also write / ~ I4(tf,tf)). In this case some explicit 
computations are possible. 

The coefficients am, A, 7 do not depend of /; the only coefficient depending on / is the first 
time of arrival x c (which is nonzero only if a < ajf). If a < clm, then 


x c = tf + a(a,p)(tf - tf) 

where a(a,p) := ’ with C ( a ’P) = f/ xp(a) ,1 d /„ p < 1/p. This means that, in the 

subcritical case, the ratio R between the length of the arrival times interval and that of the dis¬ 
turbance times interval, depends only on a and p (not on /); indeed R := (tf — x c )/(tf — tf) = 
(1 -p + C(a,p)~ 1 )~ 1 . 


The cumulative distribution function T/ can be computed using equation (3.6). Although 
when p / 0 no explicit evaluations are possible, one can see that F fl is a rescaling of the cumulative 
distribution function of the ESS obtained when the disturbance is uniformly distributed on the 
interval [0,1] (with the same parameters a and p), which we denote by F/ p . More precisely, if 
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F v = F“ ,:p is the cumulative distribution function of the ESS when / ~ then 

\tf — tj V x c ) 


(3.9) 


(recall that tf V x c = x c if a < clm while tj V x c = tf if a > aw). Using computer-aided numerical 
solutions, in Figures [5[)10| we plot the cumulative distribution functions of the ESS (solid line) and 
of the disturbance (dashed), with different parameters and disturbance intervals. In the first row, 
Figures [5j|7j we take / ~ (0.5, 0.9) while in the second row, Figures [ 8 }{Tq| we have / ~ U( 0.3,1). 

In the three figures in each row, the parameters (a,p) are, respectively, (0.2, 0.2), (5,0.2) and 
(5,0.5). This implies that the first figure represents a subcritical case (a = 0.2 > oa/(0.2) = 
3.30447). Note that in both rows in this case the arrivals start towards the end of the disturbance 
intervals and the ratio between the length of the two intervals representing the arrival times and 
the disturbance time is the same, R = 0.208552. 

The second figure of the row is a supercritical case (a = 5 > »m( 0.2)), where a fraction 
7 = 0.10142 of the population arrives at time 0. In both rows the fraction 7 is the same (7 depends 
only on a and p). 

The third figure is again a supercritical case (a = 5 > om( 0.5) = 0.941046), where strong 
competition forces a larger fraction of the population (7 = 0.456433) to arrive early. Along with 
the values of a and p, in Figures [5] |10| we write the explicit value of the maximum average fitness 
A. Note that A increases when either p decreases or a increases, the maximum competition cjm- 



Figure 8: a = 0.2, P = 0.2, a = 0.83315s. Figure 9: a = 5, p = 0.2, a = 0.120448. Figure 10: a = 5, p = 0.5, a = 0.0510315. 


As we have seen in Figures [5 10 by (3.9), it is enough to study the case / ~ 7/(0,1). In 
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Figures TTflTT) we plot the cumulative distribution function Fq ' p (solid line) corresponding to nine 
couples (p, a ), together with the cumulative distribution function of the disturbance (dashed line). 
In each row p takes the same value (p = 0.1, 0.3 and 0.5 in the first, second and third row 
respectively), while a is constant along each column (a = om( 0.1), a^f (0.3) and qm (0.5) in the 
first, second and third column respectively). In this way, figures on the diagonal represent critical 
cases, figures on the upper triangle are subcritical cases and figures on the lower triangle are 
supercritical cases. 





Figure 11 : a = clm (o.i) = 6.893, p = Figure 12 : a = clm (o. 3 ) = 2.075, p = o.i, Figure 13 : a = clm (0-5) 

0.1, A = 0.0995, 7 = 0. A = 0.3153, 7 = 0. A = 0.5123, 7 = 0. 





-c.d.f. of ESS fi 

-—— c.d.f. of disturbance 


= 0.941, p = 0.1, 


-c.d.f. of ESS /> 

-—— c.d.f. of disturbance 


Figure 14: a = a M ( o.i), p = 0.3, a = Figure 15: a = a M (o.3), p = 0.3, a = Figure 16: a = M (0.5), P = 0.3, a = 

0.0478, 7 = 0.2663. 0.296271, 7 = 0. 0.5066, 7 = 0. 





-c.d.f. of ESS // 

-—— c.d.f. of disturbance 


Figure 1 /: a — a M (o.i), p — 0.5, a = Figure 18: a — a M (o.3), p = 0.5, a = Figure 19: a — ajy(o. 5 ), p = 0.5, a = 

0.0184,7 = 0.4788. 0.2606,7 = 0.3141. 0.5011,7 = 0. 


3.2 Climate changes 

In our model the climate can be represented by the couple (p, /), that is, the distribution and the 
strength of the disturbance; more precisely the lower p the stronger the disturbance. One can argue 
that also the competition parameter a could be affected by the climate, nevertheless, in this paper, 
we prefer to think of it as a characteristic of the population. 
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When climate changes (i.e. the couple (p, f) changes), there are two interesting questions. 

1. How does the ESS change reflects a climate change? In other words, can we predict in which 
respect the new ESS will differ from the previous one? 


2. What happens if a colony keeps the same strategy of arrivals after a climate change? Does 
the average fitness of the population decrease or increase? 

Some answers can be obtained if we imagine that the climate change affects just one of the 
parameters p and /. 


The answer to question 1, is given by Theorem 3.5 If p decreases, then A (a,p) increases: after 
the population has adapted to the new climate the common fitness will have increased (of course 
supposing that the population is sufficiently large to survive the transition period). Moreover, a 
decrease of p may lead from an ESS with massive arrival at 0, to an ESS with arrivals after a date 
x c > tf. Let us discuss now what happens if / changes, that is if it moves from fa to / 2 - The fitness 
of every individual, according to an ESS adapted to fa, is A (a,p) which does not depend on fa, thus 
that will not change in case of rapid adaptation (question 1). What will change is the distribution 
of the arrival times and possibly its support. In the case of uniformly distributed disturbances, the 


answer to question 1 is given by equation (3.9): a change of the interval during which disturbances 
occur is simply reflected by a rescaling of the ESS. It is worth noting that an analogous result holds 
if we take a generic density / defined on [0,1] and we rescale it as follows 


fi(t) ■■= /( = 


t — t 


f 


f-if 


(3.10) 


the effect on the ESS in this general case is still the rescaling given by equation (3.9). 

The second question is particularly interesting in the case of rapid climate changes, since the 
adaptation of the population (moving from the old ESS to the new one) could require several 
generations; thus the changes may endanger the survival of the population. A sudden change of 
climate may be an advantage for some individuals (i.e. some arrival dates) and a disadvantage for 
other individuals. Even in the case of a simple anticipation of the disturbance the situation is 
not trivial. Suppose that fa is supported in [f ],L] and /2 is supported in [t 2 ,fa], with fa < fa. 
Individuals arriving at y > fa are sure that the disturbance is over but they might have more 
competitors still alive (for instance those arriving at y — e), thus it is not clear whether the climate 
change implies a larger or smaller fitness. The following proposition gives a partial answer: let 
us denote by 4>^\y) the fitness of an individual, in a population following the strategy p, (not 
necessarily an ESS), which chooses the arrival date y and has to face the disturbance which is 
regulated by fa. If /2 is a delay of the disturbance, that is, in the second scenario disturbances strike 
later in the season, early birds suffer a decrease of their fitness. The second part of the proposition 
tells us that, at least for uniform disturbances, if /2 delays the beginning of the disturbance but 
also anticipates the end of the disturbance season, then later birds have an advantage. 
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Proposition 3.7. Consider a distribution y and two densities f\ and / 2 on the intervals 

t*2>*2]- 


1- Ifh<t 2 then ($\y) > \y) for all y E [li,£ 2 ]. 

2. If y is the ESS associated to f\. D 2 is a random variable with density / 2 and F (t 2 < D 2 <tf) 
is sufficiently small then <j$\y) > <$\y) for all y E 

3. If f\ and / 2 are uniform densities and h < t 2 and t\ > t 2 then <f>\f\y) > 4>^\y) for all 
y>yi-= (*i ?2 - hh)/(ti +h-t 2 - ti) (where yi E [t 2 ,t 2 \)- 


Proposition |3.7| ( 1) and (3) follow from a more general result (Proposition 5.8) which allows to 
compare <f>^ and dyP in more general settings. Proposition 5.8 is fairly technical and can be found 
in Section [5} Proposition 3.7'2) implies for instance that a consistent delay of the disturbance 
(think of t 2 > t\) reduces the fitness of all individuals following the former ESS. Even if it is not 
necessary that t 2 > ti, the requirement that F(f 2 < D 2 < t±) is sufficiently small, implies that 


t'2 > ti¬ 
lt is not trivial to assess whether individuals that profit from the climate change represent a 
small or large fraction of the population following what was a former ESS. More precisely, it is 
not always clear whether the average fitness A^ increases or not. Nevertheless it is not difficult to 
provide some partial answers in the case of delay or anticipation of the disturbance. 


Proposition 3.8. Consider two densities f\ and / 2 on the intervals \t 2 ,to\ and let y be the 

ESS associated to f\. Denote by D 2 a random variable with density f 2 . 

1. a|j 2) < A^ provided that P(t 2 < D 2 <t\) is sufficiently small. 

2. a|T' ) > A^ provided that F^ftf) is sufficiently small. 


Proposition |3.8[ 1) implies that a consistent delay of the disturbance reduces the average fitness 
of the population; on the other hand, Proposition 3.8 2) implies that the average fitness of the 
population increases if the disturbance arrives so early that most of the population has not yet 
arrived when the danger is over. 

Let us discuss now the consequences of a change of p. In the previous sections we showed that 
p 1 —> X(a,p) is decreasing. This means that if p decreases and there is an instantaneous adaptation 
of the population to the new conditions (i.e. the arrival distribution y changes accordingly) then 
the fitness of each individual (hence the average fitness) increases, because more individuals will 
choose to arrive later in the season. Intuitively, if the population does not adapt to the climate 
change, a decrease of p should imply a killing of a larger fraction of the population and hence a 
lower average fitness. The following result tells us that this intuition is correct, at least for weak 
competition. 
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Proposition 3.9. Given a generic arrival strategy /a, recall the definition of the average fitness as 
Xu = X n(a,p, fjL, f) = f ^(y)/a(dy). The function p (->■ = A n(a,p, p,, f) is continuous in [0,1]. 

Moreover for all a G (0, 21og(2)] we have d p X p (a,p, p) > 0 for all p G (0,1) and the inequality is 
strict if and only if F p (x) > 0 for some x <tf. 

This means that, if the competition is low, then when the probability of surviving the distur¬ 
bance decreases then the fitness decreases as well. Roughly speaking, in case of low competition it is 
better to arrive later (and to have more competitors) than to take a chance against the disturbance. 

4 Discussion 

In this paper we developed a model for the timing of life history events when a disturbance can strike 
a population of migrators and may kill some of the individuals that incur in it. Our model therefore 
considered the biologically realistic scenario of a “soft” disturbance, i.e. an event, like extreme 
unseasonable weather, that can kill a fraction of the population, but to which some individuals 
survive. We developed our study by considering migratory birds as a model, because they are 
biological system where the timing of life history events has been intensively studied. In particular, 
we looked for the ESS for individuals that have to choose arrival time to their breeding grounds, 
and benefit from early arrival, as is often the case for migratory birds ([16|). On the other hand, 
they may incur in a catastrophe, for instance a spell of cold weather, that will kill a fraction of 
individuals that have arrived to their breeding grounds before the catastrophe. Clearly, the choice 
of focusing on arrival to the breeding grounds is purely exemplificative, and our model applies 
more generally to the timing of almost any life-history event. For example, it applies to timing of 
arrival to the wintering grounds, crossing of a geographical barrier like a mountain range, or to 
timing of reproduction, rather than arrival. From a biological point of view, the most interesting 
results we obtained can be summarized as follows. First, in presence of competition, a fraction 
of individuals arrive early during the season, in a period when they can incur in a catastrophe 
and therefore be killed. Interestingly, the stronger the competition, the earlier birds start arriving 
and may even arrive at the earliest time possible. Remarkably, there is a threshold value for the 
intensity of the competition above which a fraction of individuals arrive extremely early, i.e. at 
a time when they will certainly incur in the catastrophe. Hence, competition is able to force 
individuals to risk death if the payoff for an early arrival is sufficiently high. Second, a strong 
disturbance increases the fitness of individuals (fitness is equal for all individuals that follow the 
ESS) because, under a strong disturbance, the fraction of individuals arriving at a time when they 
can incur in the catastrophe decreases. Hence, a strong disturbance determines a later average 
arrival of the population. Actual distribution of arrival dates is therefore the balance between the 
contrasting pressures of competition and risk of death due to the catastrophe. Third, the ESS is not 
the strategy that determines the maximum average fitness of individuals in the population. Indeed, 
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a larger average fitness could be obtained if some individuals would accept a reduction in their own 
fitness for favouring other individuals. This result is not surprising since individuals are predicted 
to behave selfishly and adopt the strategy that maximizes their own fitness. We also tested our 
model under a climate change scenario, which is predicted to increase the frequency and strength of 
extreme meteorological events ([22]), like spells of unseasonal weather that can kill migratory birds. 
We obtained two further interesting results. Fourth, climate change, which is predicted to increase 
the strength of the disturbance, should determine an increase in the fitness when individuals are 
able to adapt their arrival times to the new ESS (second result above). However, the ESS implies 
a later average arrival of individuals (third result above). Differences among species or populations 
in the observed shifts in the timing of migration according to climate change may therefore be due 
to differences in the susceptibility of species or populations to extreme weather conditions, which in 
our model is accounted for by the strength of the disturbance. Since climate change is determining 
on the one side a general advancement of the timing of spring events ( [291 (33| ), but on the other 
side it also increases the frequency of cold spells in spring ( |20| ), differences in the response of bird 
species or populations to climate change should be investigated also in respect to their susceptibility 
to unseasonable weather. Fifth, if the population is unable to adjust arrival times and continues 
following the previous timing, which is no more an ESS, the fitness of individuals declines (in many 
scenarios). It has been hypothesized that many migratory bird species, particularly long distance 
migrants, may be less able than short distance ones to advance their arrival to the breeding grounds 
because they are constrained by the timing of other life-history events m>- Consequently, they are 
forced to follow an arrival strategy that differs from the current ESS, and should suffer a reduction 
in fitness. Our model therefore gives an explanation of the possible mechanisms linking response 
to climate change and population trends and explaining why bird populations that did not show 
a response to climate change are declining ([30]). There is currently debate among biologists on 
whether the observed changes in arrival dates of migratory birds can be attributed more to micro¬ 
evolutionary processes or to phenotypic plasticity ([H]). If the response to climate change is due to 
phenotypic plastic response of individuals, then probably adaptation will be fast enough to keep the 
pace of climate change. In contrast, if timing of life history events is genetically controlled, a longer 
time may be needed for the new ESS to fix in the population, and in the meanwhile individual 
fitness will be reduced. In summary, the model we developed may contribute to our understanding 
of the processes determining the timing of life history events under the biologically realistic scenario 
of a catastrophe killing only a fraction of the individuals that incur in it. Our model explains how 
competition can induce a fraction of the population to arrive very early, despite facing a higher risk 
of death, as it is documented in several species ([33]). Moreover, our model also investigated the 
effect of climate change on the timing of life history events, and demonstrated that fitness should 
decline in a scenario of increased probability of catastrophe if the population is not able to adapt 
to the new climatic conditions. 
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5 Proofs 


In this section one can find all the proofs of our results and some details about the remarks of the 
previous sections. We start with a lemma and its corollary. 

Lemma 5.1. Let {a n } n >o, {b n } n >o and {k n } n >o be such that a n ,b n > 0 for all n>0,b n >0 for 
some n > 0, k n+ 1 > k n for all n > 0 and a n = b n k n for all n > 0. Define no := min{n > 0: b n > 0} 


m : = 


En >0 a ^ n 
En>0 b n t n 


t = 0 . 


Then f is a nondecreasing function on [0, +oo). Moreover f is strictly increasing on [0, +oo) if 
and only if there exists m , n such that k rn > k n and b m ,b n > 0 . 

Proof. Note that / is continuous on [0,+oo) and differentiable on (0,+oo). We compute the 
derivative on ( 0 , +oo) as 


fit) = 


En>l 1 • En>0 En>0 ^ ' En>l nb n t h 

( En>0 b n t n \ 


En >0 Ei=o(* l)( a i+lbn—i a ibi+l) t n En >0 Ei= 0 (* + k n —i ) t n 


E^>0 b n t n 


E^>o b n t n 


0 {kn+l k 0 ){n + l) 6 n +l &0 “I - E^i =0 (* d" ^)bi+lbn—k n —i) t 


En>0 M" 


Now (k n+1 — fco)(n + l)b n +ibo > 0 and, by using j = n — i — 1, 


^ ' (i kn—i) 


n —i n —i 

+ V)bi-\-\b n —kn-i) + ^ j')bn—jbj—l(kn—j kj— l) 

2=0 j =0 

n—1 n—1 

^ 1 )^ 24 - 1 ^ 72 — 2 (^ 2+1 k n — 2 ) ^ ^ j^)bj—\b n —j{kj—\ k n —j ) 

2=0 j =0 

n—1 ^ n— 1 

^ ^ ^ 2 +i^n— 2 (^ 2 +! k n —i)(i + 1 (n i)) = — ^ ^ ^ 2+1 bn— 2(^241 k n —i)(^li + 1 71 ). 


Since (/c^+i — k n -i){2i + 1 — n) > 0 for all n,i such that n > i > 0 we have that f'(t) > 0 for all £, 
hence / is non-decreasing. 
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Moreover, if k m > k n and b m , b n > 0 then rn > n and bi+ib n -i(ki + i — k n -i){2i+\ — n) > 

bmb n (km ~ k n )(m — n) > 0 (just take i = m — 1). This implies f'(t) > 0 for all t E (0, +oo), whence 
/ is strictly increasing on [0, +oo). On the other hand if k n = k m for all n. m such that b m b n > 0, 
clearly f(t) = ko for all t E [0, +oo). □ 


Corollary 5.2. Suppose that a > a > 0 and p < 1. The the function f(t ) := ( e at — p)/(e at — p) is 
strictly increasing in [ 0 ,+oo). 


Proof. Apply Lemma 5.1 using ao = bo = 1 — p, a n = a n /n\ for all n > 1, b n = a n /n\ for all n > 1 
and k n = ( a/a ) n . □ 


The following is a brief remark which proves some of the properties of the functions 7 and A 
(the others are straightforward). 


Details on Remark \3.4\ Most of the results about the functions om, 7 , A and x c follow easily by 
checking the monotonicity and the the continuity (in each variable separately) of the l.h.s. and 
r.h.s. of the defining equations. We highlight just the main details. 


1. Since the l.h.s. of equation (3.2) is strictly decreasing, continuous with respect to p and 


strictly increasing, continuous with respect to aM and since the r.h.s. is strictly decreasing, 
continuous with respect to p we have that p <—> cim(p) is strictly decreasing and continuous. 
As for the limit we note that, for every fixed /3 E (0, + 00 ), 

r-exp(/3) A „ n if p 0 + 

if p —»• 1 “ 


P 


s; 



since l/(z 1 “ p — p) f 1/z pointwise in the first case and the Bounded Convergence Theorem 
applies, while 1 /{z 1 ^—p) t +00 pointwise in the second case and the Monotone Convergence 
Theorem applies. The limits follow easily by standard arguments. 


2 . Since the r.h.s. of equation (3.3) is strictly increasing and continuous with respect to p and 


with respect to a and since the l.h.s. is continuous and nonincreasing with respect to x c we 
have that the maps a 1 —> x c and p 1 —> x c are strictly increasing. As for the continuity, note 
that, by definition, for every (3 > x c (ao,Po, f) we have fp f f(x)dx < fj( a{) po ^ f(x)dx hence 

rh / / /- ex P( a ) 


lim ( f f(x)dx f 1 Po T f [ 

a^a Q ' J 1 3 ' ' Jl 

lim ^ f f(x)dx — (l — p + (^ j 
P^Pn J 0 J 1 


P~P0 — 


P0 


- 1 \ -1 


< 0 


exp(a 0 ) 


A~P — 


-In -1 


< 0 . 


P->P 0 Z~ r ps 

These inequalities imply that, eventually, the maximal solution x c (a,po, f) E (x c (ao,Po, /),/?) 
(resp. x c (ao,p,f) E (x c (ao,Po, /), /?))• The limits from the right can be treated analogously 
by carefully dealing with the intervals where (3 >-)• f(x)dx is constant. 
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Let us prove the limits in equation (3.5). By using the monotonicity of p H > x Cl we have that 


x c (a,p, f ) < x c (a, 0, /). Since the r.h.s. of equation (3.3) equals a/(l + o) < 1 when p = 0we 
get that x c (a,0, f) < tf. In order to compute the second limit we observe that the r.h.s. of 
equation (3.3) tends to 0 as a —> 0 + . As for the second line in equation (3.5), note that for ev¬ 
ery (3 > tf we have that fj/ f(x)dx < 1 hence fg f(x) dx— (^1 — p+ ^ ^ d /y 

fg f f(x)dx — 1 < 0 eventually as p —> pm{o) - 


3. The continuity follows from the continuity of the r.h.s. of equation (3.4) with respect to p, 
from the continuity of the l.h.s. with respect to p, a and 7 and from the monotonicity with 
respect to 7 . 

Let us prove that lim^oo j(a,p) = p. We define G(a,p , 7 ) := exp(-a 7 ) J v i d /_ p and 
we can assume p G (0,1). Clearly, ~ ife = (exp(ap) - exp(a 7 p))/p as 

a —> + 00 . Hence for every fixed 7 and p and for every e > 0 we have 

(1 + e)(exp(a(p- 7 )) - exp(a 7 (p - 1 )))/p > G(a,p,j) > (exp(a(p - 7 )) - exp(a 7 (p - 1 )))/p 
eventually as a —> + 00 . Hence 


lim G(a,p, 7 ) = 

a—>■+00 


+00 if 7 < p 
0 if 7 > p 


thus for every 7,7 satisfying 7 < p < 7 , eventually as a —> + 00 , the solution 7 (a,p) 
to equation (3.4), that is G(a,p, 7 ) — 1/p = 0, satisfies 7 (a,p) G ( 7 , 7 ). This implies 
lim a ^. 0O 7(a,p) = p. 

In order to prove that 7 (a,p) < p for all a > 0, p G (0,1), it is enough to show that 
G(a,p,p ) < 1/p for all a > 0, p G (0,1). Indeed, in that case, the solution 7 (a,p) of 
G(a,p, 7 ) — 1/p = 0 satisfies 7 ( 0 , p) G (0,p) (note that 7 >-)• G(a,p, 7 ) is strictly decreasing 
for all a > 0). Observe that 


G(a,p, 7 ) = 


exp (a(p - 7 )) - exp(—ory(l - p)) 


P 


— exp(— 07 ) 


r»exp(a) 


1 


, exp(a 7 ) ^V 1 P-p 


1 

V X ~P 


du 


exp(a(p - 7 )) - exp(-o 7 (l - p)) 

<-exp — 


exp(a) 


P 


/•ex pt 

p(-U7) / 

J exp (a- 


p 


/ , 7 ; 2(1 p ) 

P( a 7) u 


du 


exp(a(p - 7 )) - exp(— 07(1 - p)) 


p 


~ A(a,p, 7 ) 


where 


A(a,p, 7 ) := 


e _a 7 a(l — 7 )/2 if p = 1/2 

e "“ 7 2^1 (e a(2p_1) - if p ± 1 / 2 . 

When 7 = p we have G(a,p,p ) = 1/p — (e _ap ( 1_p )/p — A(a,p,p)), hence it is enough to prove 
that e~ ap< ' 1 ~ p ' ) /p — A(a,p,p) > 0, that is, F p (a) := pe ap ( 1 _p )A(a,p,p) < 1. The first case 
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is 7 = p = 1/2; thus, F p (a ) = e~ a ^ A a /8 which attains its maximum value (in [0,+oo)) at 
a = 4 and F p (2) = e~ l / 2 < 1. The second case is 7 = p/l/2 where F p (a) = p 2 (e _a ( 1_p ) 2 — 
e -ap(i-p)^ j ^ 2 p — i). Note that F p (0) = 0, F p (a) > 0 for all a > 0 and F p (a) —> 0 as a —>• +oo. 
Hence F p admits a global maximum in [0, +oo), which must be a stationary point since F p is 
differentiable. By taking the derivative F p (a) = p 2 {p{l —p)e~ ap ( 1 ~ p ' ) — (1 — p) 2 e _a ( 1 _p ) 2 )/(2p — 
1) we have that a stationary point a must satisfy g- a p( 1 _ p) — (1 — p)e~ a ^ l ~ p ^ 2 / p = 0. Hence 
Fp(a) = pe~ a ( 1 ~ p '>~ < 1. This proves that 7 (a,p) < p for all p G (0,1) and a > 0. 


Now we prove that a (-»• 7 (a,p) is strictly increasing in [clm( p ), + 00 ) for all p G (0,1). Since 
7 >-)• G(a,p, 7 ) is strictly decreasing for every fixed a > 0, it is enough to prove that a 1 - 7 - 
G(a,p, 7 ) is strictly increasing (for fixed 7 < p < 1) to obtain that a 1 - 7 - 7 (a,p) is strictly 
increasing where it is defined by equation (3.4) (by standard arguments for implicitly defined 
functions). By continuity it is enough to prove that G(a,p, 7 ) > G(a,p, 7 ) where a > a > 0 
(the case a = 0 would follow easily). 


To this aim, let us replace the variable in the integral defining G in the following way: 
z = av + (3 where a = (e“ — e“ 7 )/(e a — e“ 7 ) and f3 = —(e a+ “ 7 — e a+a 7 )/(e a — e“ 7 ) (note that 
(3 < 0 since 7 < 1). This implies that when v = e “ 7 then 2 = e “ 7 and when v = e a then 
z = e a . Whence 


P 7(a-a) 

G(a,p, 7 ) = e - “ 7 - 


= e " a7 


r” 

Je d ^ 

(*-!-(. 

J e a-y Z 1 P~p\\ 


dz 


c-^y-p-v 


z 1 p — p 
z~P\l- 


d (a-a) 




a 


dz = (*). 


Observe that (e 7 ^ a a ))/a = (e °( 1 7 ) — l)/(e a ^ 7 ) — 1). Moreover 

z'- p -p 1 - p (j) 1 - p -p/a 1 - p _ i-pO ~ (l) 1 ^ ~P(1 ~ l/a 1 ^) 

(^) l ~ p -p a ( £ ^) i -^-p a ^ c-^y-p-p 


which is strictly increasing with respect to z since z/a < (z + \{3\)/a. Hence ((z + |/3|)/a) 1_p — 
(. z/a) l ~ p is nonincreasing (since 1 — p < 1) and (z + |/3|)/a) 1_p —p is strictly increasing. This 
implies that, for all z G [e“ 7 ,e a ), 

z !-p _ „ e 7(«-0 e a(i -p) _ p e «(i- 7 ) __ 1 e “(i-p) _ 1 g a (i—7) _ 1 

— ■ - < — • __ < -. — 

^z-]3y- p a e a{l-p) _ p gd(l— 7 ) _ l — e a(l-p) _ j gfl(l-7) — 1 

e “(l-p) _ l g«(l~P) _ 1 
- e a(l-p) _ X g“(l-p) — 1 

where in the second inequality we used (x — p)/(y — p) < (x — l)/(y — 1 ) for all x > y > 1 > p, 
while in the last inequality we applied Corollary |5.2| (since 7 < p). Finally, this yields 

_ /- e “ dz 

(*) < e “ 7 / - = G(a,p, 7 ). 

J e a-y z 1 p -p 
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□ 


Remark 5.3. The continuity of the function X(a,p ) is easy. In the interval [0, om(p)] the function 
a H > X(a,p ) is strictly decreasing since the integral in the r.h.s. of equation (3.7) is strictly increasing 
for all p E (0,1). If a E [om(p),+oo) then, using equation (3.4), we have that X(a,p ) is a solution 
to 


X 


fexp(a) ^ z 


J P /\ z l ~ p ~ V 

and, since the l.h.s. of this equation is strictly increasing with respect to a and X, standard arguments 
imply that a i-)- A (a,p) is strictly decreasing. 

We show now that p i->- A (a,p) is strictly decreasing. Let us start with the first expression in 


= 1 


equation (|3.7|). In the following equation it is easy to show that the derivative with respect to p and 

exp(a) z l-p + X 


the integral with respect to z commute and 
l /-exp (a) 


dp x - 


A—p — 


p 


+ (! ~P) 


(z 1_ P — p) 2 


-dz 


fexp(a) Z 1 -p 


ln^ 1 P ) + 1 — p — (z 1 P — p) 


L (z'-P-p) 2 

r-exp(a) ^1 -p -p) + 1 _ z l-p 


> 0 


Ji {z l p — p ) 2 

for every a > 0 since the integrand is strictly positive for every z > 1; indeed the function x i-)- 
xln(x) + 1 — x = x(ln(x) — 1) + 1 is differentiable in (l,+oo) and continuous in [l,+oo) and the 
derivative is ln(x) > 0 for all x > 1. This implies that p i->- A (a,p) is strictly decreasing for every 
fixed a > 0. 

Let us consider now the second expression for A, namely pexp(-aj), which holds forp > pm(cl) 
where we recall that pm(p) is the unique solution for -:l d f_ p = ^ with respect to p. From 

equation (3.4) we have that, for every fixed a > 0, the function p i->- y(p) := pexp(—07(0, p)) is 
implicitly defined by the equation F a (y,p) = 0 where 

rexp(a) 

ip/y V l ~P~p L 

The solution to the previous equation is uniquely defined in ( Pm,P ) since y 1 —> F a (y,p ) is strictly 
increasing for every fixed p E (0,1) and F a (pM,P ) < Pm/pm — 1 = 0, F a (p,p) = p/pm — 1 > 0. We 
can compute, using similar arguments as before, 

1 f*M a ) v i-p \ n (v) + l 


F a (y,p) = V 


d P F a (y,p ) = - 


(p/yf-P -p +y Jp/y 


> - 


1 


+ 


1 


{v l ~P — p) 2 

pexp(a) ^ 


0 p/y) 1 p -p 1 ~p 

1 1 

+ 


y 


(p/y ) 1 p -p 1 -p 


J p/y 
> 0 


X-p - 


-dr 


-dr 


P 
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where in the last equality we used the fact that, since F a (y,p) = 0, then y f p Jy^ V i-p_ p ^ v = 1- The 
last inequality holds since p > y. By standard arguments, d p F a (y,p ) > 0 implies that p *->• A (a,p) 
is strictly decreasing for every a > 0 . 


Before proving Theorem 3.5 we prove that an ESS is of the form p 
is an absolutely continuous probability measure. 


7 $o + (1 — 7 )v, where v 


Lemma 5.4. Let p he an ESS and fix p E [0,1). Then for some 7 £ [0,1], p = 7<5o + (1 — 7 )v, 
where v is an absolutely continuous probability measure. Moreover [minsupp(/z),i/] D supp(/i) D 
[minsupp(/i),fy] HEsupp(/) and z/min supp (/it) < minsupp(i') then supp(i/) D Esupp(/). Finally, 
0 P is constant on [minsupp(//), + 00 ) and equals A. 


Proof. We already noted that we can write equivalently 


r»+00 


(j) fl (x) = exp(~aF p (x)) / exp(a(l - p)F p (z))f(z)dz + p / f(z)dz 


since Esupp(/) C [0, tf]. From this equation we can see easily that 0 is right-continuous on [0, + 00 ); 
moreover it is left-continuous at x if and only if p({x}) = 0. If (a, 0] is such that p((a,0]) = 0 
then (a,0] n Esupp(/) / 0, if and only if 0(0) > 0(a). Indeed, 


0n(0) ~ 0/j,( a ) = exp(-aE /i (a))^exp(a(l -p)F fl (a)) - p'j j f(z)dz>0. (5.11) 

As a consequence, if p((a, 0)) = 0 then (a, 0) n Esupp(/) / 0 if and only if lim t _^- 0y,(t) > 0 p (a). 
Moreover if {x n } ne N is such that x n E supp(^), x n / x for all n E N and x n fx then x E supp(/it) 
and p({x}) = 0. The first assertion, namely x E supp(^), comes from the fact that the support of 
a measure is a closed set. Suppose, by contradiction, that x is an atom; clearly 0(x n ) = A for all 
n, since p is an ESS, and lirn^^.- 0 p (t) > 0 fl (x) hence 0 p (x) < A which is a contradiction. This 
proves that the only atom, if any, of p must be 0. Thus p = jdo + (1 — 7)17 where 7 = p( 0) and v 
is nonatomic. 

We prove now that if a < 0, a E supp(p) and 0 E Esupp(/), then 0 E supp(/r); this implies 
supp (p) D [minsupp(/it), tf] n Esupp (/). By contradiction if 0 0 supp (p) there exists e > 0 such 
that (0 — e, 0 + e) n supp( / u) = 0. Let 0 := max (supp (/it) n [0, 0\) then 0 < 0. Since 0 E Esupp (/), 


fp-ii f( x )dx > 0. Clearly (0,0 + e) n supp(/j) = 0, as a consequence of equation (5.11), 0 < 


lim f _^ +£ )- 0 p (t) — 0n(0) = lim i _ 5 .( / 3 +£ )- 0 p (f) — A. Hence there exists t E (0 — e, 0 + e) such that 
0fi(t) > A but this contradicts the definition of ESS. 

Let us prove that [min supp (p), tf] D supp (/it). Observe that, for every 0 > tf we have 


0^(0) - 0n(tf) = (exp (-aF p (0)) - exp (~aF p (t f )) / exp(a(l - p)F p (z))f(z)dz < 0. 

■Jo 
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On the other hand, if /3 £ supp(ft), there exists e > 0 such that F^ffl + e) > F^(tf) and f3 + 
e £ supp(ft). This implies + s) < 4>/i(tf) which contradicts the definition of an ESS. Thus 
[minsupp(ft),t/] D supp(ft). 

If rninsupp(/i) < minsupp(t/) then 7 > 0 and 0 = minsupp(ft) < min supp(i'). On the other 
hand, supp(z/) C supp(ft) C supp(z^) U {0}, that is, supp(ft) \ {0} = supp(^) \ {0}. Since when 
minsupp(ft) = 0 we have supp(ft) D Esupp(/), then supp(^) \ {0} D Esupp(/) \ {0} which implies 
supp(i/) = supp(t/) \ {0} D Esupp(/) \ {0} = Esupp(/). 


We are left to prove that v is absolutely continuous. First of all we note that, from equation (2.1), 
= exp (-aF^(x))H^(x) where 

rx /»+oo 

H fl (x)= exp(a(l-p)F fl (z))f(z)dz + p f(z)dz 
■JO Jx 

is absolutely continuous and nondecreasing on [0, if]. Since the only atom of ft, if any, is 0 and 7^ is 
right-continuous, clearly (f>n is continuous on [0 ,tf]. It is well-known that any open set in M can be 
decomposed into a disjoint, at most countable union of open intervals, hence (0, +00) = (supp(ft) \ 
{0}) U U jejlj where {Ij}j £ j is an at most countable disjoint union of open intervals in (0, +00). 
By definition, <^(x) = A for all x £ supp (ft). Since supp(7/) D (minsupp(ft) n Esupp(/),i/) then 
j T f{x)dx = 0 for all Ij C (minsupp(ft), +00). As a consequence of equation ( 5.1 1[ ) , (j)^ is constant 
on such intervals Ij, but, since the value of q ^ at the extremal points of Ij is A, then 0„(x) = A for 
all x £ Ij. This proves that <f>n(x) = A for all x > minsupp (ft). Hence A = exp(— aF fl (x))H fJj (x) for 
all x > minsupp(ft) which implies 


F n(x) = < 


-log 

a 


H u( x ) 


if x > min supp (ft) 
if x < min supp (ft). 


(5.12) 


By composition, F) is clearly absolutely continuous on (minsupp(ft),ff). Hence F v is absolutely 
continuous on (minsupp(ft),i/) (since F^{x) — (1 — 7 )F u (x) = 7 for all x > 0) and this implies that 
v is an absolutely continuous measure (since supp(t/) C [minsupp(ft), t/]). □ 

Given a measurable subset I C M, we denote by T x (/) the set of real, measurable functions on 
I which are integrable with respect to the Lebesgue measure. When a different measure p needs to 
be specified, we write L 1 (J, p) instead. 


Remark 5.5. A well-known characterization which will be used in the sequel is the following 
Lebesgue’s fundamental theorem of calculus (see for instance f3E\ . Theorem 7.18]). A function G is 
absolutely continuous on a compact interval I if and only if there exists a function g £ T 1 (/) such 
that for some for all) a £ I and for all f3 £ I, G(/3 ) — G(a) = J'f g{x)dx; in that case, for 

almost every x £ I, G is differentiable at x and G'(x) = g(x). 
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This implies that an absolutely continuous function G on I is constant if and only if G'{x) = 0 
for almost every x E I. 

Moreover consider two functions h : I —» J, k : J —>• Y; if k is absolutely continuous and h 
is monotone and absolutely continuous then k o h is absolutely continuous. Similarly, if k is a 
Lipschitz function and h is absolutely continuous then ko h is absolutely continuous. If h\ : I —> M 
is absolutely continuous then h + h± is absolutely continuous and the same holds for h ■ h\ if I is 
compact. 

Finally ifk is differentiable everywhere and h is differentiable almost everywhere clearly we have 
(k o h)'(x) = k\h(x)) ■ h'(x) almost everywhere and, if in addition k oh is absolutely continuous, 


k o h{x) — k o h(xo) = / (koh)'(z)dz = / k'(h(z)) ■ h'(z)dz = (*) 


(5.13) 


’x 0 


’x 0 


for all xq,x 6 I. If, in addition, h is absolutely continuous then 


rHx) 
'h(x 0 ) 


k\y)dy = k o h(x) — k o h(x o) = (*). 


(5.14) 


A function k is locally Lipschitz on I if and only if for every x$ £ I there exists 5 > 0 and M > 0 
such that \x — xo| <5, \x — xo\ < 6 implies \h(x) — h{x)\ < M\x — x\. By elementary analysis, if h 
is a locally Lipschitz function on I and I' C / is compact then there exists M > 0 such that for all 
x, x E I' we have \h(x) — h(x)\ < M\x — x\ (that is, h is globally Lipschitz on every compact subset 
of I). It is easy to show that a locally Lipschit function on I is absolutely continuous on I (one can 
prove it on every compact subset and then use the fact that I is the union of an increasing family 
of compact subintervals). Hence the following result holds. 

Given an interval I, if k is locally Lipschitz on I, differentiable everywhere and h is absolutely 


continuous on I then the equalities (5.13) and (5.14) hold; moreover k o h is constant on I if and 
only if (kf o h)(x) ■ h'(x) = 0 almost everywhere in I. 


Proof of Theorem \3.S\ By Lemma 5.4 we know that an ESS can be written as y = 7 < 5 o + (1 — 7 )^, 
where v is an absolutely continuous measure and supp(/r) D [x^Hf] n Esupp(/), where x^ := 
minsupp(/x) (we will show later that x^ = x c ). Since F^(y) — (1 — 7 )F u (y) = 7 for all y > 0, we 
have that F( is absolutely continuous on [ 0 ,+ 00 ). 

Since fn is absolutely continuous, we can take the derivative of fn in (2.1), for almost every 

V > 0 , 

^(y) = -a(l - 7 )y(y)(fu(y) + exp (-aF^y)) exp(a(l - p)F li (y)) - p f(y), (5.15) 

where g is the derivative of v. From now on it will be tacitly understood that the derivatives and 


equalities involving them, are defined and hold almost everywhere. From Lemma 5.4 we know that 
for all y > Xfj, we have 4>fi(y) = A hence (f'Ay) = 0. Thus 

f(y ) 


g(y) = 


o(l -y)A 


exp {-aF^y)) exp(a(l - p)F fl (y)) - p 
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Let z := exp (aF^). Using the last result of Remark 5.5 (by taking I = [0,+oo), k(x) := exp (ax) 
and h := F) t ), the previous equation is equivalent to 

z\y)= fi f(z(y) 1 - p -p), 

which, since z(y) l ~ p > p for all y > 0 , is in turn equivalent to 

*(y) = f(y) 

z(y) 1 ~ p — p A 

this is equivalent to 


with the condition z(x fl ) = exp (ay). By Remark 

/•exp(a(l—'y)F„(x)+a'y) 


5.5 


d^ 


yl~P — 


P 


f f(y)dy, 
J Xu 


(5.16) 


J exp(a7) 

which, once we prove that x c = x^, is equivalent to equation (3.6). This proves that given 7 , F v is 
uniquely defined by the previous equation. From the previous equation, for all x G [x^,tf], 

fexp(a(l—7)i r I /(x+e)+o7) ^ ^ r x + £ 


L 


^_ rQ 

exp (a (1—7)^(2—e)+<27) % ~ P 


1 rx- t-e 

- j f(y)dy, 


which implies that supp(^) = [x^,,tf] n Esupp(/). Now, since supp(/r) \ {0} = supp(^) \ {0}, from 
Lemma 5.4 we have that either p = v (hence supp(/r) = [x^, tf] n Esupp(/) and minsupp(/x) = 
minsupp(zv)) or 7 > 0 and x^ = 0. In this case supp(^) = Esupp(/) and supp(/r) = Esupp(/)U{0}. 

In particular, since x^ G supp(/j), when p = v (that is, 7 = 0) we have that x M G Esupp(/). 
Clearly the measure v is supported in [0, tf] , but it is a probability measure if and only if F v (tf) = 1, 
therefore, by using equation (5.16), x M and 7 must satisfy 

/*exp(a) 


dz 


/exp(a7) Z P p 


1 /•*/ 
A L 


f{y)dy. 


(5.17) 


We note that A = </> A 1 (minsupp(/^)) = (p^(x /Jr ). If 7 = 0 then p = v hence, from equation (5.12), 
by continuity of F^ at x = minsupp(ju) = x^ we have A = that is 


A = 


r x n 

/ f{y)dy 

Jo 


+ P [ f(y)dy = 1 - (1 -p) f f(y)dy 


(5.18) 


since exp(a(l — p)F^(y)) = 1 for all y G [0, x^). 


If 7 > 0 then x M = 0 and, from equation ( 2 . 1 ), A = <^(0) = pexp(— 07 ) which is the second 
line of equation (3.7). 


dz 


r*exp(a) 


Let a < clm- If 7 > 0 then by equation (|5.17|) and the discussion thereafter, 

l rexp(a M ) 

p 

> 


h z 

/•exp (a) 


1 -v — 


> 


p 


'1 


Z 


dz 

1 -p — 


p 


1 1 

= V > 


I exp(a7) Z p p A p 
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which is a contradiction. If 7 = 0 then by plugging the explicit value of A = 1 — (1 — p) f* f f(y)dy 
into (5.17), we get (3.3) (with instead of x c ) and since we showed that G Esupp(/), thus it 
must be the (unique) maximal solution to equation (3.3), that is x^ = x c . This proves that a < clm 
implies 7 = 0 and x c > tj. Since the r.h.s. of equation (3.3) is strictly less than 1 (resp. equal to 
1) when a < clm (resp. a = om) then equation (3.3) yields x c > tj, (resp. x c = tf). In order to 
obtain the first line of equation ( |3.7[ ) just consider the expression of f* f f(y)dy given by the r.h.s. of 
equation (3.3) and plug it in equation (5.18) (recalling that = x c ). 

Let a > clm- If 7 = 0 then, by the definition of ajvf, equation (5.17) and the discussion thereafter, 

Sl[ f(y) d v 


1 r*f 

A L 


f(y)dy = 


1 - (1 - p) fJ f{y) d y p Jl 


/*exp (a) 


< 


r */ 

dz 

yl~P - 


1 

< - = 


rexp(a M ) 


?1—P — 


P 


P 


1 r*f 

a L 


f(y)dy, 


which is a contradiction. If 7 > 0 then x^ = 0 (which coincides with the definition of x c given 
before equation (3.3)) and by plugging the explicit value of A = pexp(— 07 ) in equation (5.17) we 
have 

pexp(a) /-exp(a(l— 7 )) 


- = exp(— aj) 


L 


V " " 7exp(a 7 ) ^ V - p J 1 (zexp(aj)) 1 P-p' 

which has a unique solution 7 . Note that the r.h.s. of the equation is a decreasing continuous 


function of 7 , say iL(y); moreover it(0) > 1/p by (3.2) and 

/■exp (a M ) 


K(1 - a M /a) = 


dz 


rexp(a M ) 


(zexp(a - aw )) 1 p -p 


< 


dz 


yl~P — 


P p 


This implies that the unique solution 7 £ (0,1 — aM/a)', 7 < p was proved in Remark 5.3 


□ 


Before giving the details on Remark 3.6 we need a Lemma and another remark. 

Lemma 5.6. Let p and v be two finite measures on M. Then the following conditions are equivalent: 

(1) f h(x)p(dx) > f h{x)v{dx) for every nondecreasing, measurable h G L 1 (M,/r) nL 1 (M, is); 

(2) p((x, + 00 )) > v((x, Too)) for every iGK and p(M) = i/(M); 

(3) f h(x)p(dx) < f h(x)u(dx) for every nonincreasing, measurable h G L 1 (M, p) n L 1 (M, u); 


(4) p{{~ 00 , x]) < i/((— 00 , x]) for every iGK and p(R) = i/(M). 

If one of these conditions holds we write p^v. 

Moreover the inequality in (1) is strict if and only if there exists y such that /r(/i _1 ([y, Too))) < 
^(/i _1 ([?/, Too))) (that is, if and only if there exists y such that p(/i _1 ((y, Too))) < v (h" 1 ((y,+ 00 ))) ). 

Finally the inequality in (3) is strict if and only if there exists y such that p(/i _1 ([?/, Too))) < 
zv(/i _1 ([?/, + 00 ))) (that is, if and only if there exists y such that p(/i _1 ((y, Too))) < ;/(/i~ 1 ((//, Too)))/. 
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Proof. The equivalence between (1) and (2) (or between (3) and (4)) is a classical result of measure 
theory: it is a slight modification of the arguments in [431 Section l.A.l]. Moreover, to prove that 
n(R) = i'(M) just take condition (1) (or condition (3)) and consider first h = 1 and then h = —1. 
Under /r(M) = z'(M) the equivalence between conditions (2) and (4) is trivial. In particular the 
equivalence between the previous four conditions hold even if we take the set of nonnegative (or 
nonpositive) measurable functions instead of L 1 (M, y) n L 1 (M, v). 

Finally note that y H > y{[y, Too)) and y H > v{\y, +oo)) are continuous from the left while y H > 
y((y, Too)) and y e-» i/((y, +oo)) are continuous from the right; this implies that there exists y such 
that /r(/i _1 ([?/, Too))) < v(h~ l (\y, Too))) if and only if there exists y' such that /i (h —1 ((?/, Too))) < 


Remark 5.7. Given a random variable X with law y, then the composition £ := F fl o X between the 
variable and its cumulative distribution function is a random variable with values in Rg(iT) C [0,1] 
such that P(£ < t) < t for all t E [0,1]; moreover P(£ < t) = t if and only if t E Rg(F fl ) := 
F m (M). More generally P(£ < t) = sup {F^r): F^(r) < t}, thus if y(z) > 0 then P(£ < t) = 
lim s ^,- F^s) = F^z) - y{z) for all t € [lim s ^ 2 - F^s), F^(z)) = [F^z) - y(z), F^z)). 

To be precise if we consider {x: y{x) > 0} =: {xf. i E J}, the at-most-countable set of discon¬ 
tinuity points of F^ (where J C NJ, and define I := (J iej[F/j.(xi) ~ //(xj), F^xf)) we have 


W) 


t if t fL I 

Fn{xi) - n{xi) ifte [i^(xj) - n(xi), Ffj{xif). 


This implies in particular that if F^ is continuous then the random variable £ is uniformly distributed 
on [0,1] otherwise, according to Lemma \5. 4 £ >~ Unif(0,1) (that is, the law of ( stochastically 
dominates a the uniform distribution on [0, 1\). 


Details on Remark \3.(A For a generic /, by performing a time rescaling, we can assume without 
loss of generality that if = 1. Consider the following family of absolutely continuous (with respect 
to the Lebesgue measure) measures / u n (■) = n\ ■ n[l — 1/n, 1]| with support in [0,1]. (where | • | 
represents the Lebesgue measure on M). If we define c n (f,p ) := / 0 ' f(z)dz + p fi_ 1 / n f(z)dz 
then c n (f,p) t 1 as n —> oo and we have 






= —on(y—l+l/n) 


(y)lin(dy) = j 

> C e~ an( ' y ~ 1+1 / n ' ) c n (/, p)ndy = c n (f,p) [' e~ ay 'dy' f 1 ~ CXp( ~ a 

J1—1/n Jo 


Cn(f, P ) T / ^-PMz-l+l/n) - p) f( z)dz \ ndy 
Jl—l/n 


l— 1 /n JO a 

as n —> oo (according to the Monotone Convergence Theorem). This example can be modified in 
many ways: for instance, one can take a family of absolutely continuous measures, with strictly 
positive densities on [ 0 , 1 ], which approximate the previous densities nl[in the L 1 norm (we 
do not give more details on this). 
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We are left to prove that j cf) ^(y) y(dy) < (1 — exp (—a))/a for all y. We start by proving that 
for every measurable, nonincreasing function h and every z G M we have +oo ^ h(F fl (y)) y(dy) < 
f(F (z) 1 ] (and, if h is strictly decreasing on [0,1], then the inequality is strict if and only if F ^ 


is not continuous on [z, +oo)). By stochastic domination (see Remark 5.7), if we denote by P^ the 
law of £, f R h(Fn(y))n(dy) = Jj () ^ h(s)P^(ds) < Jj 0 ^ h(s)ds and, for a strictly decreasing function, 


the inequality is strict if and only if T), is not continuous on M (apply Lemma 5.6 to P^ and to the 
uniform distribution on [0,1] and choose y := h(s) for some s G [lim„,^ 2 - F^(w), F^{z)) whenever 
this interval is nonempty; note that, for a strictly decreasing function, h~ 1 ([h(s),+oo)) = (—oo,s] 
for all s ). Take zgR; then 




/i(s)P|(ds) = 




K F ^{y))^y) = / h(Ff,{y))y(dy) (5.19) 

J (£,+00) 


since 


F n 1 ((F fl (z), 1]) D (z, Too) and y(F^ 1 ((F fl (z), 1]) \ (z, +oo)) = 0. 

By the previous equation we just need to prove that f^ F ^ ^ /i(s)P^(ds) < f^ F ^ ^ h(s)ds and 
this can be done in two different ways. 

) ,Ffj,(z)]) equals F^(z), that is, the 


The quick way is to note that, according to Remark 5.7 


lebesgue measure of the interval [0,, F^(z)]; hence the measure P^ restricted to (F^z), 1] dominates 


the Lebesgue measure on the same interval, thus Lemma 5.6 applies. 

An alternate proof is as follows; if we define 

%) := ('* W if S - 

\h(F F (z)) if s<F f (z) 

then h is nonincreasing hence, by stochastic domination, f R h(F fl (y))/j,(dy) < f, () ^ h(s)ds. Thus 


I (jZ,+Oo) 


h(Ff,(y))y{dy) = 


K F ^y))^y) = / H F ^(y))^y) - / K F n{y))v(dy) 

J] R J (—00,2:] 


f (z,+oo) 

= [ K F u.{y))^y) - H F n(z))n((-oo,z]) 

Jr 

< [ h(s)ds-h(F,j,(z))iJ,((-oo,z])= [ h(s)ds - h(F /J ,(z))F IJ ,(z) 

J[ 0,1] J[ 0,1] 

= f h(s)ds — f h(s)ds 

•7 [0,1] J[ 0,F„(*)] 

= / h(y)&y = / Hy)dy 

J(F^(z), 1 ] J(F^z), 1 ] 


(5.20) 

and, if h is strictly decreasing on [0,1], the inequality is strict if and only if F ^ is not continuous 


on (z, + 00 ) (since h is strictly decreasing on [F^z), 1] we use equation (5.19) and we can apply 
again Lemma |5.6| to Pg and to the uniform distribution on [0,1] by choosing y := h(s ) for some 
s G [lim, ( .^ 2 - F^(w), F^(z)) whenever this interval is nonempty and 7 G (z + 00 )). 
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Clearly 


rX 


rt f 


<My)Mdy) = E exp(-a£) / exp(a(l - p)E„( 2 :))/(z)dz + p / /(z)dz 


L JO 


/x 


Hence, using Fubini’s Theorem and the inequality (5.20) with h(s) := e as , we have 

rv 


0 M (y)/x(dy) = e aF ^ y) 


I o 


(e “( 1 F f v( z ) — p)f(z)dz + p fi(dy) 


= p f e aF ^ p{dy) + f (e ^ 1 P ^ F>1 ^ — p)f(z) f e aF ^p(dy)dz 

J J «/ (z,-|-oo) 

p [ e~ ay dy+ [(e a ^~ p)F ^ z) -p)f(z) [ e~ ay dydz 

J fo.n J J(f„(z), ii 

+ J{e a ^- p)F ^ z) -p)f{z)- 


'[ 0 , 1 ] 

1 -e" a 


v-aF^z) _ -a 


= P 

= P 


a 

l-e~ a 

a 


+ J(e~ apF ^ -pe~ aF ^)f{z) 


1 _ e -a(l-F M (z)) 


dz = (**) 


and the inequality is strict if and only if F ^ is not continuous. Now (1 — e~ a ^ t '^ z ^)/a < (1 — e~ a )/a 
and e~ apFfi ^ — pe~ aF/1 ^ < 1 —p (since s i —y e~ aps —pe~ as is decreasing in [ 0 , +oo)); moreover each 
of the previous inequalities become a strict inequality if and only if F^z) > 0. Whence 


, . 1 - e~ a , , 1 - e~ a 

(**) < p - b (1 -p) - 

a a 

and there is a strict inequality if and only if inf supp(/z) < supEsupp(/). □ 

Proof of Proposition \ 3. 7[ 1. Applying Proposition |5.8[ we have easily that f y fi(x)dx > Jf f 2 (x)dx 

for all 2 < y, y <E [t 1 ? t 2 ]> hence 4$\y) > 4$\y)- 

2. If t 2 > t\ the conclusion follows from Proposition |3.7| (1) . Let us suppose t 2 <ti- Now let us 
evaluate (pf (y) for y E [t 2 ,t i]: 

^\y ) = exp(-aF M (y)) f exp(a(l - p)F f _ l (x))f 2 (x)dx + pexp(-aF f ,(y))F{D 2 G [y,t i)) 

“ t.2 

+ pexp(-a)P(-D 2 > ti) = (*). 


Note that this equality holds even if t± > t 2 since in that case ¥(D 2 > ti) = 0. 

(*) < exp(-apF fl (y))P(D 2 E [t 2 ,y]) + pexp(-aF lx (y))F(D 2 € [y,t i]) + pexp(-a)P(L > 2 > h) 
< P(L »2 G [t 2 ,i i]) +pexp(-a)P(L > 2 > ?i) 

which is smaller than exp(— a) if P(-D 2 G [t 2 ; i l ]) is sufficiently small. When p is an ESS, then 
for all y E [i[, t\] we have (j$\y) = A(a,p) > exp(—a). If t 2 < t ] then the proof is complete. 
If t 2 > f l then the conclusion follows using Proposition |3.7|(1) for y E [t\ , f 2 ] • 
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3. 


Consider the solution y\ to the equation (y — — t\ ) = (y — t 2 )/{t 2 

{t 2 ti—t 1 t 2 )/(t 2 +ti — ti — t 2 )', clearly for all y > y\ we have (y — tfij/fti — ti) 


thus Proposition 5.8 implies 4>^\y) > <j$\y)- 


~t 2 ), namely y 1 : = 


□ 


Proposition 5.8. Consider two probability densities f\ and f 2 . Let us define 4>^ andfiffl according 
to equation (2.1) using f\ and f 2 respectively. Fix y > 0 such that J^ fi(x)dx > 0 or f Q v f 2 (x)dx > 

0 . 


1. Define the set of strong maxima from the left of as M^(y) := {z E [0, y]: F/j,(s) < 
F^(z) for all s < zj. If 


j fi(x)dx > J f 2 (x)dx, Mz E M^(y) 


(5.21) 


then 4$\y) > <t$\y)- In this case (f^\y) = <f>\f\y) if and only if equality holds in equa¬ 
tion (5. 2l| ) for all z E M /i (y). 

2. If fy fi(x)dx > Jq f 2 (x)dx and a is sufficiently small then 4>^\y) > (f\f\y). 


Proof of Proposition\5.8[ According to equation (2.1) we have 


^(y) = exp(-. aF^ (y))[ [ (exp(a(l - p)F tl (x)) - p)ffix)dx + p 

o 


, * = 1 , 2 . 


1. We just need to prove that J^ (exp(a(l -p)F li (x)) -p)fi{x)dx > / 0 y (exp(a(l - p)F fi (x)) - 
p)f 2 (x)dx. This follows easily from the fact that x >-)• exp(a(l— p)F fl {x))— p is a nondecreasing 
function and applying Fubini’s Theorem. Indeed, for every z > 0, the set {x > 0: exp(a(l — 
p)F fl (x)) — p > z} is an interval of type [xo(z),y\ (for some xo(z) E [0, y]); clearly xq( z) = 0 
for all z E [0,1 — p\. Hence 


(exp(a(l -p)F fl (x)) -p)ffix)dx = f \ f 

Jo l J 0 

-I 


y r r ex p( a P~p)Fn(x)) —p 


d z 


fi(x) dx 


rexp(a(l—p)Ffj,(y))—p . 
3 

rexp(a(l—p)F fl (y))—p r 


L J{xG R: exp(a(l—p)F Al (x))— p>z} 

dz 


fi(x)dx 


dz 


j 

J\Xi 


(5.22) 


= (!-£>) / fi{x)dx + 
A 0 , 3 /] 


fi(x)dx 

L J[x 0 (z),y] 

exp(a(l—p)F M (2/))—p 


1 ~P 


i 


[ 3 : 0 ( 2 ), y] 


fi{x) dx 


dz. 


The equivalence of the equalities is as follows. The “if” part is trivial. As for the re¬ 
verse implication, note that 2 ^ xq (z) is left-continuous and that s = xq (z) for some 
z E [0,exp(a(l — p)F fl (y)) — p] if and only if s = inf{t > 0: Fyft) > a} for some a < F^(y), 
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that is, if and only if xq(z) 6 M M (y). If fp fi(x)dx > fp f 2 (x)dx then $\y) > <$\y) 


follows from equation (5.22). If J'P(, r) fi(x)dx > f 2 (x)dx for some xq(z) > 0 then by 

the continuity of w *->• jp f\ (x)dx and the left continuity of z >-)• xq(z) there exists e > 0 such 
that fi(x)dx > fz 0 ( s ) f' 2 (x)dx for all s £ (z — e, z]\ again equation ( |5.22 ) yields the strict 
inequality $\y) > $\y)- 


2. It follows from equation (5.22) 


$Hy) ~ <f>jj?Hv) > (1 “ P) [ (/i(®) - f 2 (x))dx - (exp(o(l - p)Ff,(y )) - 1). 

*'[ 0 , 2 /] 


by taking the limit as a goes to 0. 


□ 


Proof of Proposition 3.8. 1. It follows from the first part of Proposition |3.7| 


2. If t\ > t 2 then, by Remark 3.6, = (1 — exp(—a))/a, which is larger than Aj^ and the 

proof is complete. If t x < t 2 , write A^ = (1 — exp(—o))/a — e. Note that, if y > t 2 , then 
(p'n (y) = exp(— aF^{y)) f t 2 exp(a(l — p)F fl (x))f 2 (x)dx. Thus, if we integrate with respect to 

/i: 


> [ f>P\y)dp{y) 

Ju 


= ~(exp(—aF(t 2 ) — exp(—a )) 


r *2 


f 2 {x)dx+ / exp(a(l - p)F tl (x))f 2 (x)dx 


' t-i 


> -(exp (~aFJt 2 )) - exp(-a)) > -(1 - exp(-a)) - £ 
a a 

where the last inequality holds if Fpift 2 ) is sufficiently small. 


□ 


Proof of Proposition^ 79| By using some basic results of mathematical analysis, it is easy to check 


that d p f 4> fl (y)y(dy) = f dpf p(y)p{dy) where, from equation (2.1), 


ry rtf 

dp(j>n(y) = exp(-aF M (y)) -a exp(a(l - p)F lx (x))F tl (x)f(x)dx + / f(x)dx 

JO Jy 


From equation (5.20) we have 


[ e~ aF ^ y) fi{dy) < [ e~ ay dy = 

J(z, + 00 ) J(F„.(z), 11 


D -aFp,(z) _ „—a 


1 ] 


(5.23) 
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According to Fubini-Tonelli’s Theorem 


dpX p — dp 


= —a 


4>ii(y)^(dy) = / 9 p ^(y)/i(dy) 


+ 


J tf e aa-p)F»(x) Fix{x) f {x )(J 

, ,,/w (X 


e -aF„ 


e -“^ 


(i,+oo) 

dx 


dx 


/(*) 


' [o,®] 


e - a ^ (y V(dy) - ae a ^- p)F ^ x) F^x) 


e~ aF » 


(x,+oo) 


dx 


Let us study the expression between brackets by means of the function h p defined by the following 
equation 


e ~ a F 


'(x,+oo) 


h p (x) := f e- aF ^n(dy)-ae a ^- p)F ^Fp{x)( [ 

J[0,x] W(; 

> [ e~ aF ^ y) n(dy) - e~ apF ^ x) Fp(x)( 1 - e -«(i-M*))) =: h p (x) 

J [0,a:l 


which holds for every x G [0 ,tf], p G [0,1] (the inequality comes from equation (5.23)). Clearly 
p i->- h p (x) is nondecreasing hence h p (x) > ho(x) for all x G [0,ty]. In particular p H > h p (x) 
is strictly increasing if F p (x) > 0; in this case h p (x) > ho(x). Suppose that a < 21og(2); thus 
e -az _|_ e -a(i-z) > 2e _a / 2 > 1. Now, since F p { 0) = //(0), 

h 0 (0) = p(0)e~ a ^ - fi(0)(l - e -«(!-M°))) = ^(o) ( e _a M0) + e -a(i-A*(o)) _ i) > 0 . 
Moreover, if z > x > 0 we have 

h 0 (z) - h 0 (x) = [ e~ aF ^p(dy) - Fp(z)( 1 - e -«(i-^*(*))) + F„(x)( 1 - 
= [ e- aF ^p(dy)-(Fp(z)-Fp(x))(l-e- a ^- F ^) 

J(x,z\ 

+ F p (x)(e~ a ^- F ^ z)) - e-^i-^d))) 

> (F p (z) - F p (x))(e~ aF ^ + _ i) 

+ F p {x){e~ a{ - 1 - F ^ z)) - e - a d-^d))) > o 

where the last inequality is strict if and only if F p (z) > F fl (x) (we also used the fact that 
f(x z] e ~ aFlx ^y-(dy) > e~ aFft ^p,((x, z]) = e~ aF ^ z \F p {z) — Fp{x))). This means that, for all 
x G [0,t/], p G [0,1] we have h p (x) > h p (x) > h o(x) > ho(0) > 0. Observe that d p X p(a,p, y, f) = 
fo f f(x)h p (x)dx. Whence, for all a < 2 log(2), for all and for all /, we have that d p X tl (a. p, /j, /) > 
0 for all p G [0,1]. 

Observe that F p (x) = 0 for all x < tf if and only if tf > inf supp(/i). Hence, if F fl (x) = 0 for 
all x <tf clearly 4>p{y) = 1 for all y G supp(/u) which implies A p(a,p, y, f) = 1 for all p G [0,1]. 
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On the other hand if F p {x) > 0 for some x < tf we have that h p (z) > 0 for every z 6 [x,tf], 
p G (0,1]; thus for all a < 21og(2), for all p and for all /, we have that d p \ p (a,p, p, /) > 0 for all 
PG( 0,1], 

□ 


References 

[1] L. Belhadji, D. Bertacchi, F. Zucca, A self-regulating and patch subdivided population, 
Adv. Appl. Probab. 42 n.3 (2010), 899-912. 

[2] D. Bertacchi, N. Lanchier, F. Zucca, Contact and voter processes on the infinite percolation 
cluster as models of host-symbiont interactions, Ann. Appl. Probab. 21 n. 4 (2011), 1215-1252. 

[3] D. Bertacchi, G. Posta, F. Zucca, Ecological equilibrium for restrained random walks, 
Ann. Appl. Probab. 17 n. 4 (2007), 1117-1137. 

[4] D. Bertacchi, F. Zucca, Critical behaviors and critical values of branching random walks on 
multigraphs, J. Appl. Probab. 45 (2008), 481-497. 

[5] D. Bertacchi, F. Zucca, Characterization of the critical values of branching random walks on 
weighted graphs through infinite-type branching processes, J. Stat. Phys. 134 n. 1 (2009), 
53-65. 

[6] D. Bertacchi, F. Zucca, Approximating critical parameters of branching random walks, 
J. Appl. Probab. 46 (2009), 463-478. 

[7] D. Bertacchi, F. Zucca, Recent results on branching random walks, Statistical Mechanics and 
Random Walks: Principles, Processes and Applications, Nova Science Publishers (2012), 289- 
340. 

[8] D. Bertacchi, F. Zucca, Branching random walks and multi-type contact-processes on the 
percolation cluster of Z rf , Ann. Appl. Probab. 25 n. 4 (2015), 

[9] K. Bessho, Y. Iwasa, Variability in the evolutionarily stable seasonal timing of germination 
and maturation of annuals and the mode of competition, J. Theoret. Biol., 304 (2012), 66-80. 

[10] C. Both, A.V. Artemyev, B. Blaauw, R.J. Cowie, A.J. Dekhuijzen, T. Eeva, A. Ene- 
mar, L. Gustafsson, E.V. Ivankina, A. Jarvinen, N.B. Metcalfe, N.E.I. Nyholm, J. Potti, 
P.A. Ravussin, J.J Sanz, B. Silverin, F.M. Slater, L.V. Sokolov, J. Torok, W. Winkel, J. Wright, 
H. Zang, M.E. Visser Large-scale geographical variation confirms that climate change causes 
birds to lay earlier. Proceedings of the Royal Society of London B 271 (2004), 1657-1662. 

[11] J.M. Calabrese, W.F. Fagan, Lost in time, lonely, and single: Reproductive asynchrony and 
the allee effect, American Naturalist, 164 (2004), n.l, 25-37. 

[12] N. Champagnat, R. Ferrire, S. Meleard, From individual stochastic processes to macroscopic 
models in adaptive evolution, Stochastic Models 24 (2008), n.l, 2-44. 


36 



[13] B. Chan, R. Durrett, N. Lanchier, Coexistence for a multitype contact process with seasons, 
Ann. Appl. Probab. 19 (2009), 1921-1943. 

[14] A. Charmentier, P. Gienapp, 2013. Climate change and timing of avian breeding and migration: 
evolutionary versus plastic changes, Evolutionary Applications, 7 (2013), 15-28. 

[15] P.A. Cotton, Avian migration phenology and global climate change, Proceedings of the Na¬ 
tional Academy of Science USA, 100 (2003), 12219-12222. 

[16] H.Q.P. Crick, D.W. Gibbons, R.D. Magrath, Seasonal changes in clutch size in British birds. 
Journal of Animal Ecology, 62 (1993), 263-273. 

[17] R. Durrett, S. Levin, The Importance of Being Discrete (and Spatial), Theor. Pop. Biol., 46 
(1994), n.3, 363-394. 

[18] W.J. Ewens, Mathematical population genetics. I. Theoretical introduction. Second edition. 
Interdisciplinary Applied Mathematics, 27. Springer-Verlag, New York, 2004. 

[19] O. Garet, R. Marchand, R. B. Schinazi, Bacterial persistence: a winning strategy?, Markov 
Process. Related Fields , 18 n. 4 (2012), 639-650. 

[20] L. Gu, P.J. Hanson, W.M. Post, D.P. Kaiser, B. Yang, R. Nernani, S.G. Pallardy, T. Meyers, 
The 2007 eastern US spring freeze: increased cold damage in a warming world?, BioScience, 
58 (2008), 253-262. 

[21] IPCC 2012, Managing the Risks of Extreme Events and Disasters to Advance Climate Change 
Adaptation, A Special Report of Working Groups I and II of the Intergovernmental Panel 
on Climate Change [C.B. Field, V. Barros, T.F. Stocker, D. Qin, D.J. Dokken, K.L. Ebi, 
M.D. Mastrandrea, K.J. Mach, G.K. Plattner, S.K. Allen, M. Tignor, and P.M. Midgley 
(eds.)]. Cambridge University Press, Cambridge, UK, and New York, NY, USA, 582 pp. 

[22] IPCC 2013, Climate Change 2013. Cambridge University Press, Cambridge. 

[23] Y. Iwasa, S. A. Levin, The timing of life history events, J. Theor. Biol., 172 n. 1 (1995), 33-42. 

[24] Y. Kang, N. Lanchier, The role of space in the exploitation of resources, Bull. Math. Biol. 74 
(2012), 1-44. 

[25] N. Lanchier, The role of dispersal in interacting patches subject to an Allee effect, 
Adv. Appl. Probab. 45 (2013), 1182-1197. 

[26] N. Lanchier, C. Neuhauser, Stochastic spatial models of host-pathogen and host-mutualist 
interactions I, Ann. Appl. Probab. 16 (2006), 448-474. 

[27] N. Lanchier, Contact process with destruction of cubes and hyperplanes: Forest fires versus 
tornadoes, J. Appl. Probab., 48 (2011), n.2 352-365. 

[28] S. Meleard, Some stochastic models of interacting diffusion processes and the associated prop¬ 
agation of chaos, Proceedings, Stochastic Modelling in Biology, Heidelberg (1990), 107-125, 
World Scientific. 

[29] A. Menzel, T.H. Sparks, et al. European phenological response to climate change matches the 
warming pattern, Global Change Biology, 12 (2006), 1969-1976. 


37 



[30] A.P. M0ller, D. Rubolini, E. Lehikoinen, Populations of migratory bird species that did not 
show a phenological response to climate change are declining, Proceedings of the National 
Academy of Sciences of the USA 105 (2008), 16195-16200. 

[31] A.P. Mpller, W. Fiedler, P. Berthold, Effects of climate change on birds, Oxford (2010), Oxford 
University Press. 

[32] S. Muller, Recurrence for branching Markov chains, Electron. Commun. Probab. 13 (2008), 
576-605. 

[33] I. Newton, Weather-related mass-mortality events in migrants, Ibis 149 (2007), 453-467. 

[34] I. Newton, The ecology of bird migration, (2008) Academic Press, London. 

[35] R. Pemantle, A.M. Stacey, The branching random walk and contact process on Galton-Watson 
and nonhomogeneous trees, Ann. Prob. 29, (2001), n.4, 1563-1590. 

[36] E. Post. S.A. Levin, Y. Iwasa, N.C. Stenseth, Reproductive asynchrony increases with envi¬ 
ronmental disturbance, Ecology, 55 (2001) n.4, 830-834. 

[37] T.A. Richter, The effects of density-dependent offspring mortality on the synchrony of repro¬ 
duction , Evolutionary Ecology, 13 (1999), n.2, 167-172. 

[38] W. Rudin, Real and complex analysis, Third edition, (1987), McGraw-Hill, New York. 

[39] A. Satake, A. Sasaki, Y. Iwasa, Variable timing of reproduction in unpredictable environments: 
Adaption of flood plain plants, Theoret. Pop. Biol.,60 (2001), n.l, 1-15. 

[40] N. Saino, R. Ambrosini, Climatic connectivity between Africa and Europe may serve as a basis 
for phenotypic adjustment of migration schedules of trans-Saharan migratory birds, Global 
Change Biology 14 (2008), 250-263. 

[41] N. Saino, R. Ambrosini, et al, Climate warming, ecological mismatch at arrival and population 
decline in migratory birds, Proc. R. Soc. B (2011) 278, 835-842. 

[42] R. Schinazi, Mass extinctions: an alternative to the Allee effect, Ann. Appl. Probab., 15 
(2005), n. IB, 984-991. 

[43] M. Shaked, J. G. Shanthikumar, Stochastic orders, Springer Series in Statistics (2007), 
Springer, New York. 

[44] M.D. Schwartz, R. Alias, A. Aasa, Onset of spring starting earlier across the Northern Hemi¬ 
sphere, Global Change Biology, 12 (2006), 343-351. 

[45] A.M. Stacey, Branching random walks on quasi-transitive graphs, Combin. Probab. Corn- 
put. 12, (2003), n.3 345-358. 

[46] G. Walther, E. Post, P. Convey, A. Menzel, C. Parmesan, T.J.C. Beebee, J. Fromentin, 
O. Hoegh-Guldberg, F. Bairlein, Ecological responses to recent climate change, Nature 416 
(2002), 389-395. 

[47] F. Zucca, Survival, extinction and approximation of discrete-time branching random walks, 
J. Stat. Phys., 142 n.4 (2011), 726-753. 

[48] F. Zucca, Persistent and susceptible bacteria with individual deaths, J. Theor. Biol., 343 


38 



(2014), 69-78. 


39 



